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Abstract 

Top quark pairs are produced with strongly correlated spins in the partonic reactions 
qq ti and gg ti. A complete description of these effects in terms of the spin density 
matrix of the tt system in leading order QCD is given. We further discuss the prospects 
to observe the spin-spin correlations at pp and pp colliders by measuring suitable angular 
correlations among the t and t decay products. 



1 Introduction 



After the recent discovery of the top quark |^ , the detailed study of the properties of this 
particle will be a major subject of experiments at the (upgraded) Tevatron and at future 
colliders. An intriguing feature of the top quark is that due to its heaviness it decays on 
average before it can form hadronic bound states. Moreover, the typical spin flip time is 
much larger than the lifetime of the top . Thus a possible polarization of the tt system 
induced by the production mechanism will be transferred to its decay products. Due to 
the dominant parity violating decay mode t — > Wb the t and t "self-analyze" their spins. 
The spin information may be extracted by forming angular correlations among the t and 
i decay products, thus allowing for a variety of tests of the standard model (SM) and 
extensions thereof (see, e.g. ||^- [p^). 

While at hadron colliders the longitudinal polarization of top quarks is practically 
zero due to parity invariance of QCD, a nonzero polarization transverse to the production 
plane is induced by absorptive parts at the one-loop level. The prospects to observe this 
order ag effect at the Tevatron and the Large Hadron Collider (LHC) have been studied 
in detail in |21| (see also Apart from this single quark polarization, the t and 

i are produced with strongly correlated spins in qq ti and gg ti, which are the 
dominant partonic production processes at the Tevatron and the LHC, respectively. In 
fact, these spin-spin correlations are of order one at the level of the partonic reactions. It 
is the aim of this letter to discuss the prospects to unravel these effects. An experimental 
verification of the feasibility to extract spin information on top quarks would clearly be 
important for any further proposals to study quantities related to the top spin. 

In the next section we will give a general description of the ti spin state in terms of 
a spin density matrix. QCD induced spin-spin correlations may be described in general 
by four "structure functions". The leading order results for these functions are given. 
Observables built from the spin operators of t and t allow to discuss the magnitude of the 
correlations at parton level. In Section 3 we will construct observables which are directly 
measurable (on an event by event basis) in pp, pp tiX with subsequent ti decays. 
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2 The spin density matrices for qq tt and gg tt 

Although the spin of an unstable particle produced in high-energy reactions is no directly 
observable quantity, it is very useful to introduce the concept of a spin density matrix for 
the tt system. Using the narrow width approximation for the top quark, we may view the 
reactions considered here as the production and subsequent decay of on-shell top quark 
pairs. The spin information may then be extracted on a statistical basis from the decay 
products of the t and t. 

We first discuss the reaction q{pi) + q{p2) t{ki) + t(/c2), where the momenta refer 
to the partonic cm. system, Pi + P2 = 0. The complete spin information is encoded in 
the (unnormalized) spin density matrix W^, 

= E m,a,)tik2,f3i)\r\qip,),qip2)r 

C colors, spins 

{t{h, a2)t{k2, /32)|r|g(pi), q{p2)). (1) 

Here, a, (3 are spin indices, A'",^ denotes the number of colors, p = Pi, k = ki and the 
sum runs over the colors of all quarks and over the spins of q and q. The factor 1/4 ■ 
takes care of the averaging over spins and colors in the initial state. The matrix structure 
of W in the spin spaces of t and i is 

= A«ll ® U + ■ (T ® U + ■ U ® cr + Cf^a' ® , (2) 

where a* are the Pauli matrices and the first (second) factor in the tensor products refers 
to the t (t) spin space. (The spin operators of t and t are simply given by (Tjl® 11 
and 11 ® cr/2.) For a detailed discussion of B? and its symmetry properties, see p!6| . 
We will now specialize on reactions mediated by strong interactions. Imposing P and CP 
invariance on Rq and using rotational invariance, we are left with the following structures: 

= 



(3 dz 

p X k 

\p X k\ 

Cij = cl{s,z)6ij + cl{s,z)piPj + cl{s,z)kikj + cl{s,z){pikj + (4) 



Bt = Bl = bl{s,z)^-^ (3) 



Here, /3 = yl — Amf/s is the velocity of the top quark in the partonic cm. system, 
p = p/\p\, k = k/\k\, z = p ■ k is the scattering angle and s = {pi + ^2)^ is the partonic 
cm. energy squared. (The notation, in particular the numbering, is adopted from [p!6| .) 
For gg —>■ ti we have an analogous decomposition with l/A^c* ~^ '^/{^c ~ 1)^ (S)- The 



functions b 



3' 



q,g derive from absorptive parts in the scattering amplitude pTl. Here 



we are interested in the spin-spin correlation functions Cq, c\, c\ and c\. They are already 



induced at Born level. The leading order results can be found in [|16[. We give them here 
in a compact form for completeness. 



qq tt 

2 2^r-l 
Denne Kq = tt — — . 

A" 



Co 



64 



C5 



Then 



K,{2 + {z'-l)(3') 



2Kq(3^ 1 



z^P^ 



-2k 



zf3' 



99 tt 

Using the abbreviation Kg 
= 2k 



a 
2k 



J) 
C4 



„9 
C5 



l + 2p\l-z^){l-(3^)-p^z^ 

2\ n2 



we have 



AKg{l - z')P 



z{l - z^)(3^ 



-Ak, 



(5) 



(6) 



Due to Bose symmetry, R^{p, k) = R^{—p, k), which means that C0 4 5 are even func- 
tions of z whereas cq is odd in z. In leading order, the structure functions for qq —>■ ti 



have the same symmetry properties, because the reaction proceeds through a single virtual 
intermediate gluon, implying that also is invariant under p —p. 

We may now construct a list of "observables" built from the spin operators of t and t: 

d\ = <T ® a, 

di = pa0p(T, 

Oi = k - a ®k - (T, 

di = {p-CT®k-(T + p-(T®k-(T)/2 {i = q,g). (7) 

At the level of tt production from partons the expectation values of these quantities for 
fixed s and z are defined as 

where the trace is over the spin spaces of t and i. They are given by linear combinations 
of the functions Co(s, 2), . . . , Cg(s, z), divided by the unpolarized differential cross section 
A^{s,z). In particular, we find 

m = 1, (9) 

which can be easily understood from the fact that the ti pair is produced from a single 
spin one boson in qq collisions at Born level. At threshold, one can show that the quantum 
numbers of the top quark pair are given by for qq tt and ^Sq for gg — > tt . 
In particular, 

lim(af) = -3, (10) 

which can be verified using (^. As a further example we plot in fig.l the rapidity distri- 
bution of (Os), defined by 

j\dztT{R^dl)5{n-T;i ^^^^ 



with 



j\ dz ii{Ri 
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at a value j3 = 0.5. We choose the rapidity of the top quark in the partonic cm. system 
as variable here instead of the scattering angle for later comparison with observables for 
pp, pp tlx . . . , for which rapidity is a convenient variable. 

The expectation value {0\) corresponds to a helicity correlation studied by Stelzer 
and Willenbrock ||2^ and also by Mahlon and Parke ||22|, namely 



_ da\tRtR + tLtL)/dz - dd-^tjjtL + tLtji)/dz 
~ da\tRtR + tLiL)/dz + da\tRiL + tLiR)/dz' 

Here, the indices L and R correspond to left- and right handed particles, respectively, 
and da^ denotes the partonic cross section, i = q,g. We find agreement with the values 
given in for this correlation integrated over the scattering angle and folded with the 
parton distribution functions both for Tevatron and LHC energies (+40% and —31%, 
respectively). Rather than to specialize on finding ways to trace this correlation of he- 
licities, we will in the following try to extract as much information on the spin density 
matrix as possible, i.e. construct angular correlations among the ti decay products in 
close correspondence to all four "partonic" observables defined in (^. 

3 Angular correlations for semileptonic tt decays 

The spin-spin correlations discussed in the previous section must be traced in the decay 
products of the t and i. The spin information in the production density matrices is 
transferred to the decay products through the parity violating decays of the top quarks. 
For semileptonic decays t i'^Ui+b, the normalized decay spin density matrix p of the 
top quark in leading order (and using the narrow width approximation for the W boson) 
reads (see, e.g., |T^) 

f+ 0+ h\ Qxf+{1-Xf+) dxi+dQi+ 
Pit - t.,.b) = [11 + q.. ■ -] (14) 

where u = m'^/m^, Xi+ = 2E^+/mt G [u, 1] is the scaled energy of the lepton, and c{l+ is 
the direction of the lepton. The asterisk refers to the rest system of the decaying quark. 
For hadronic decays, t — > W~^b — > bqq', we have, if we use the tagged b quark as spin 
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analyzer, 

pit W+h ^ hqq') ^'^ ^'^ ~ ^ 



11 H q! ■ cr 



(15) 



47r 

The corresponding decay spin density matrices p for the top antiquark are derived from 
the above ones by the replacements xi-, q|+ "Q^- and dQi+ dfli- in ([l^) 

and by — > — q^, cif^fe — >■ (iilj in (|l^). We may also use the W boson as spin analyzer 
instead of the b quark. We get the corresponding decay density matrix from (^) by 
ql —Qw+y dQy[r+. In particular, the W boson has the same spin analyzer 

quality as the b quark. The expectation value of any observable constructed from the 
momenta of the final state particles will involve a trace over the spin spaces of t and i 
of the form tr K^pit X) (g> pit X') . Thus the spin information of the production 
density matrix i?* is recovered in the decay products in a statistical sense. 

We will concentrate on decay channels where either the t or the t decays leptonically 
and the other quark decays hadronically, i.e. on the decay modes 

t W^b bqq', 

i rue-b (16) 

and the charge conjugated ones. These decay modes turn out to be especially suited to 
construct observables which are sensitive to spin-spin correlations: The charged lepton is 
the most efficient spin analyzer (cf. QT^) ), while in the same event the momentum of the 
top quark (or antiquark) may be reconstructed from its hadronic decay products. For 
nonleptonic decays, since charm tagging is difficult, we will use the bottom quark as spin 
analyzer. (Alternatively, we may use the W boson.) This leads to a suppression factor 
(1 — 2uj)/{l + 2uj) ^ 0.43 in all quantities we will discuss. One may also consider double 
leptonic decays where this suppression factor is absent. However, apart from forming 
only ~ 1/9 of all tt decays, these events do not allow for a reconstruction of the t and/or 
t rest system on an event-by-event basis due to the unseen neutrinos. Any correlation 
constructed from the laboratory momenta of the two charged leptons suffers from a large 
"background" from the unpolarized cross section, i.e. of eqs. (^, |]). In contrast, the 
quantities we will use get contributions from spin-spin correlations only. 



In close correspondence to the observables (0) we define, for tlie decay modes ([T6|): 



01 

04 



{ql-Pp)iqe- -Pp), 
{qI ■ kt){qe- ■ kt), 

ql ■ Pp){qi- ■ kt) + {ql ■ ktMe- ■ Pp)] /2. 



(17) 



Here, quantities witliout an asterisk are defined in the laboratory frame, carets denote unit 
vectors and pp is the beam direction. The top quark momentum kt has to be reconstructed 
from its hadronic decay products for a measurement of Oi, . . . O4. The momentum of the b 
quark was boosted into the top quark rest system]^. Analogous observables (5i^2,3,4 may be 
defined for the charge conjugated decay modes. The observable O4 gives zero if integrated 
over a symmetric rapidity interval in our leading order calculation because, as discussed 
above, R^'''{p, k) = i?^'^(— p, k). We therefore also define 

a = sign(n)04, (18) 



with 



1 ( Et + Pp-k; 

r+ = - in 

2 



(19) 



We evaluated the correlations {Oj) (j = 1,2,3,5) for pp collisions between 1.6 and 4 

TeV and for pp collisions between 8 and 16 TeV with nit = 180 GeV. We found only a 

weak dependence on the choice of the parton distribution functions. In the results below, 

we used the parametrization with Q'^ = Ami. applied cuts on the top quark 

transverse momentum \k^\ and rapidity: For the lower energies, we used \k^\ > 15 GeV, 

\rt\ < 2, for the higher energies, the cuts \kj\ > 20 GeV, \rt\ < 3 were imposed. 

*A note of caution: The top quark rest frame defined in equation and ( [l5| ) is different from the 
one defined in (0); the former is defined through a rotation-free boost from the partonic cm. system 
(where the matrix i?* is defined), while the latter is related to the hadronic cm. system through a pure 
boost. Thus they differ by a Wigner rotation which has to be taken into account in the theoretical 
calculation of the expectation values of the observables ( p7| ) . 
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We will first discuss the case of pp collisions. The results are shown in fig. 2. The 
correlations are largest for small cm. energies; there the correlation {O2) reaches a value 
~ 3.5%. The analogous correlations for the charge conjugated decays of the top quark 
pair have exactly the same values due to CP invariance. It might seem surprising that 
the effects are quite small remembering that we had spin-spin correlations of order one at 
the level of ti production from partons, cf. d). The suppression comes about as follows: 
As mentioned before, we lose a factor of 0.43 by using the b quark (or W boson) as 
spin analyzer. Moreover, integrating over the directions of the b quark and of the charged 
lepton generates roughly a factor of 1/9; thus the magnitude of an angular correlation built 
from these directions is ~ 5%. This is the price we have to pay for using observables which 
can be measured on an event-by-event basis and are strictly zero in the absence of spin-spin 
correlations. We will now consider signal-to-noise ratios in order to estimate the statistical 
significance of the correlations. The statistical fluctuations of our observables are given 
by AOj = ^{O]) - {Oj^. We find AOi ^ 0.58 for all energies (since (Cf) = 1/3) and, 
at = 1.8 TeV, AO2 ^ 0.36, AO3 ^ 0.36, AO5 ^ 0.25. Measuring the analogous 
correlations for the charge conjugated decays increases the statistical sensitivity. For 
example, if Nhi- denotes the number of 6-tagged, reconstructed events of type (|16D , and 
we have the same number of events in the charge conjugated channel, we get a statistical 
significance S2 for the combined correlations {O2) and {O2) of 

S.-^^g2^\/^«0.14V^ (20) 

at y/s = 1.8 TeV. In order to establish the spin-spin correlation at the 3cr level, we 
would therefore need Nhi- ~ 500, which is in reach of the upgraded Tevatron. For the 
other three correlations shown in fig. 2, we find the following statistical sensitivities at 
y/s = 1.8 TeV, again combining the decay modes (|TB|) with the charge conjugated ones: 
^1 ^ 0.076v^, ^3 ^ O.OTTVA^, ^5 ~ O-IS^K^- 

In fig. 3 we show the rapidity distributions {Oi6{rt — r'^) , . . . , {O^dirt — r[)) for 
a/s = 1.8 TeV. Note the similarity of {0^5{rt — r^)) (dotted line in fig. 3) and the 
corresponding rapidity distribution ([TI|) of the correlation {0\) in fig. 1. This similarity 
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in form is to be expected, since the partonic process qq ti dominates at Tevatron 
energies. It also holds for the other correlations and their counterparts in qq ti. 

We now turn to pp collisions. All four correlations in this case only depend weakly on 
the pp cm. energy. Between ^/s = 8 — 16 TeV they take the values {Oi) ~ (—2.37) — 
(-2.23)%, {O2) ^ (-0.13)-(-0.24)%, (C3) ~ (-0.72)-(-0.68)%, and (C5) ~ (-0.25)- 
(—0.33)%. Note especially the smallness of {O2), which in the case of pp coUisons around 
^/s ~ 2 TeV was the most sensitive correlation. We again get the same numbers for 
the charge conjugated decay modes. At = 14 TeV, we also determined the statistical 
fluctuations of the observables; as mentioned before, AOi ~ 0.58, and the other values are 
AO2 ~ 0.42, AO3 ~ 0.38, AO5 ~ 0.29. Assuming again an equal number of events in the 



decay channel (16) and the charge conjugated one, we get by combining both correlations 



the following statistical sensitivities at the LHC: Si ~ 0.055\/ N^i-, S2 ~ 0.007^/Nhi^, 
Ss ~ 0.025^/Nl,i- , and S5 ~ 0.016^y N^- . Here it is useful to consider also 

O, ^ (p, X ql) . {p^ X q,_) = - O2. (21) 

Since = (O^) — (Og)) for the combined correlations {Oq) and {Oq) a statistical 



sensitivity of Sq ~ 0.073VA^f,^- at the LHC. Note that in Oq effects of the boost from the 
partonic cm. system to the laboratory frame drop out. 

The effects are smaller at the LHC than at the Tevatron, but we have many more 
events. For example, assuming Nf^i- = 10^ = Ni^+, we can establish the spin-spin corre- 
lations to 7.3cr by measuring Oq and Oq. In fig. 4 we show the rapidity distributions of 
the correlations at ^/s = 14 TeV. Due to the dominance of gluon fusion at LHC energies, 
the shapes of the curves are similar to the corresponding rapidity distributions of the 
correlations at parton level (Ci,2,3,4) defined in (0). 

4 Conclusions 

We have shown that the sizable spin-spin correlations induced at leading order QCD 
in qq ti and gg — > ti may be measured both at the Tevatron and the LHC. For such 
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measurements, angular correlations among the decay products in "semihadronic" ti decays 
are especially suited. The correlations we propose get nonzero contributions only from 
tt spin-spin correlations — which are predicted by perturbative QCD. All calculations 
were performed at leading order and using the narrow width approximation for the top 
quark. At present, NLO corrections are only known for the functions A'^'^ in (13), 



which determine the production rate of tt pairs |25|. For a more refined theoretical study 
of the spin-spin correlations, we would need the complete spin density matrix in NLO as 
well as an estimate of the effects of non-factorizable contributions p6[. In summary, the 



experimental study of spin-spin correlations of top quark pairs seems very interesting and 
feasible. 
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Figure Captions 

Fig. 1 Rapidity distribution (|ll]) of {Ol) ai (3 = 0.5. 

Fig. 2 Correlations (Oi) (full line), (C2) (dashed line), (O3) (dotted line), and (C5) 
(dash-dotted line) as a function of the cm. energy for pp collisions. 

Fig. 3 Rapidity distributions {Oi5{rt — r[)) (full line), {02S{r-t — r[)) (dashed line), 
((D'i5{rt — r[)) (dotted line), and {046{rt — r[)) (dash-dotted line), at ^/s = 1.8 
TeV for pp collisions. 

Fig. 4 Same as Fig. 3, but for pp collisions at ^/s = 14 TeV. 
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